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^ , We present a large scale survey of life-times for spontaneous fission in the regime of 

Ch ■ super-heavy elements (SHE), i.e. nuclei with Z=104-122. This is done on the basis of 

the Skyrme-Hartree-Fock model. The axially symmetric fission path is computed using 
a quadrupolc constraint. Self-consistent cranking is used for the collective masses and 
associated quantum corrections. The actual tunneling probability is estimated by the 
WKB approximation. Three typical Skyrme forces are used to explore the sensitivity of 
the results. Benchmarks in the regime Z=104-108 show an acceptable agreement. The 
general systematics reflects nicely the islands of shell stabilization and the crossover from 
o-decay to fission for the decay chains from the region of Z/N=118/176.. 



1. Introduction 

Self-consistent mean-field (SCMF) models have proven to be re liable tools for de- 
scribing nuclear structure and dynamics, for reviews see l ^ l ^ l They are steadily 
developing to improve the quality of the predictions and to accommodate more and 
more observables. One of the very demanding observables are fission life-times. Fis- 
sion barriers have been discussed already in the early stages of SCMF models and 
even been used as benchmark for calibration SI. Fission life-times are much more 
involved as their computation requires not only the potential energy surface along 
the fission path, but also the corresponding collective masses and a safe estimate of 
the collective ground state correlations for the init ial s tate. There are thus not so 
many self-consistent calculations of fission life-timei^^ - mostly using still approxi- 
mate masses and quantum corrections (=zcro point energies^). The vast majority 
of calculations employs the microscopic-macroscopic method which c ombi nes shell 
corrections with a macroscopic liquid-drop model background, see e.g. l^iSI [ s the 
aim of this contribution to explore the systematics of fission life-times all over the 
landscape of super-heavy elements. This is done using the Skyrme-Hartree-Fock 
method as one widely used nuclear SCMF model 

The paper is outlined as follows: In section [2] we briefly review the theoretical 
background and explain how we compute fission life-times. Results are presented 
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in section [3] covering tests with experimentally known life-times, comparison with 
a-decay life-times and a large scale systematics. 

2. Formal framework 

The starting point for the self-consistent microscopic description is the SHF energy 
functional 

E = £ , SkyrmoO,T, J;j,CT,x;x) (1) 

which is expressed in terms of a few local densities and currents obtained as sums 
over single-particle wave functions: density p, kinetic density t, spin-orbit density 
J, current j, spin density er, and pair density x where each occurs twice, once for 
protons and once for neutrons. The zero-range pairing functional is handled in a 
stabilized BCS ansatz to achieve smooth transitions of the occupation numbers. For 
details see^l There exist various parameterizations for the Skyrme functional. In 
order to explore the possible sensitivity of the fission life-times to the parameteriza- 
tion, we confine the survey to three sufficientl y di fferent parameterizations: SkP as 
a force with effective nucleon mass m */m = lQ31 SkI3 as a fit which has very low 
mass m* /m = 0.6 and maps the relativistic iso- vector structure of the spin-orbit 
force ^21, and Sly6 as a fit which has m*/m = 0.7 and includes information on 
isotopic trends and neutron matter ^1. 

The computation of the fission p ath an d the ingredients of the corresponding 
collective Hamiltonian is detailed in | 14 | 15 | gj ve a ]-, r i e f summary. The mean- 
field equations are derived variationally from the given energy functional. They are 
complemented by a quadrupole constraint to generate the fission path thus reading 

[£-A&o]|*,)=£|# g ) , q=($ q \Q 20 \® q ) (2a) 

V(q) = E SkyImc (p q ,T q , Jq,Xq) , (2b) 

where h is the mean-field Hamiltonian and p q is the local density for the state 
similarly for r q , J g , and \q- The optimal fis sion p ath should, in fact, be generated 
by adiabatic time-dependent Hartree-Fock LlLLhI The quadrupole constraint is a 
plausible and generally used approximation to that. To compute the collective mass 
along the fission path, we need to explore the dynamical response of the system 
to changing deformation, commonly called self-consistent cranking. To that end 
we use the collective-momentum operator P20 (and its collective momentum p) as 
additional constraint. The P20 is deduced as generator of deformation, the response 
to P20 creates a momentum-dependent path and subsequently total energy, from 
which the mass is finally obtained as second order term in collective momentum, 
altogether 

Ao: P2o\$ q ) ocid q \$ q ) , (3a) 
[h - XQ20 - mAo] |*gp) = £\®q P ) , (3b) 
B2o = ld 2 P2( E( Pqp ,...)\ . (3c) 
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Analogously, the momentum of inertia Q x for rotations about the x- and y-axis is 
computed by self-consistent cranking. The potential V as given in eq. (|2b|) has to be 
augmented by quantum corrections to account for angular-momentum projection 
and spurious zero-point quadrupole motion ^ yielding the net potential 

V ^ v-V-B^-^^-^L . (4, 

Having determined the properly corrected collective potential V(pt2o) and in- 
verse collective mass B(q), we compute the fission life-time in the standard semi- 
classical fashion using the WKB expressions for tunneling probability W and rep- 
etition time T: 

W^expf-2/ dqJ V ty~ E ) , T = H [ dq—. 1 . (5) 

V Jb V S 20 / J a y/B 20 (E-V{q)) 

The integrals are evaluated by trapezoidal rule with separate handling at the di- 
vergent \f E — V at the classical turning points. A quantity which enters very crit- 
ically these expressions is the tunneling energy E, which is also the ground state 
energy of the nucleus before fission. That is computed with great care in a fully 
quantum-mechanical treatment of the collective quadrupole oscillations using all 
five q uadrupole degrees-of-freedom, for details of that part of the calculations see 



nve q v. 



It is to be noted that the fission path employs only axially symmetric shapes. 
One knows from actinides that triaxial shapes can produce barriers which are 0- 
2 MeV lower. The present results are thus to be understood as an upper limit 
on barriers and lifetimes. Actinides show the famous double-humped barrier. The 
outer barrier dissappears for SHE leaving only one, the formerly inner barrier 
This barrier is related to still symmetric shapes while asymmetry develops quickly 
when stepping down further towards fission. The fact that we have only one barrier 
simplifies the survey of fission properties of SHE. 



3. Results and discussion 

Fig. [T] shows collective potential and mass along the fission path of 260 Rf . The 
quantum corrections on the potential (lower panel) reduce the barrier by an im- 
portant amount making 3 to 6 orders lower life-times than would be obtained with 
the raw potential One recognizes a steep increase of the energy at deformation 
a 2 o ~ 0.6 connected with a quick change of hexadecapole momentum on®. The 
system moves from one ridge (with high «4o) to another ridge (with low a^o ). The 
seg regation of the fission landscape into ridges is well known for actinides, see e.g. 

still visible for that smaller SHE, and much smoothened with further increasing 
system size. The quadrupole mass (upper panel) shows large fluctuations. High val- 
ues are related to regions of level crossings. The detailed pattern are quantitatively 
important and can hardly be simulated by some constant collective mass. 
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Fig. 1. Example of poten- 
tial energy (lower) and col- 
lective mass (upper) along 
the axially symmetric fis- 
sion path of 260 Rf . For the 
potential energy, we show 
the raw expectation value 
V and the corrected one V, 
sec eq. Q- 
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Fig. 2. Fission life- 
times computed with 
three different Skyrmc 
paramctcrizations, as 
indicated, and com- 
pared with data from 
I21I22I23I24I25I26I 
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As a benchmark, wc compare in fig. [2] with experimental results available at the 
lower edge of SHE. There are large differences in the predictions from the various 
forces. But the experimental values are nicely bracketed by the band of predictions 
and the trends within each force compare well with experiment. The deviations of 
±4 orders of magnitude look huge. But they seem bearable in view of the extremely 
delicate balance in computing the tunneling rates. 

We are now going to present systematics of results all over the landscape of 
SHE where we consider all elements which are found to be stable against immedi- 
ate nucleon emission at ground state and along the whole fission path. Fig. |3] shows 
the fission bar riers. They agree in value and systematics with the previous survey 
of barriers^! One sees two islands of enhanced fission stability (high barriers), 
one of deformed SHE around Z/N=104/152 and the other of spherical SHE around 
120/184. The magic neutron number N=184 is also clearly marked while no unam- 
biguous sign of a magic proton number can be found in that upper island ^1 All 
three forces agree in predicting a ridge of very low barriers between the upper and 
the lower stability island. Barriers however, are mainly indications of stability. The 
full information on path, mass, and quantum corrections is required to calculate 
life-times. 
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Fig. 3. Fission barriers for a broad variety of super-heavy elements and for three Skyrme param- 
eterizations as indicated. 



Fig.[4]shows the corresponding fission life-times. The difference in barrier heights 
from to 12 MeV translates to life-times from almost immediate decay to 10 16 s, 
demonstrating again the enormous sensitivity of fission life-times to any ingredient 
in its computation. The long lived SHE are found in the two islands of stability. 
Practically immediate decay appears in the ridge of instability between the islands. 
The variance of the predictions is moderate in the islands (see also fig. [2]) and 
significantly larger in the unstable region. The latter region is not only unstable 
but also very hard to control theoretically. 

The competing decay channel for many SHE is a-decay. We have compu ted 
the a life-times from the Q a reaction energies using a recently improved recipe ^Sl 
based on the Viola systematics. Fig. [5] shows the results. The a life-times change 
much more smoothly. There are no shell effects visible and the variance of life-times 
all over the landscape is very much smaller. Comparing with fig. |4j one sees that 
a-decay is indeed the dominating decay channel in the islands of stability. Stepping 
down from the upper island, one sees a crossover of a-decay to spontaneous fission 
at abo ut Z =112 or 110, much in accordance with the experimental findings so far, 
see e.g. H^l. 

The relation between Q a v alue and a life-time is comparatively simple when 
adopting the Viola systematics It is interesting to ponder whether one could 
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Fig. 4. Fission life-times for a broad variety of super-heavy elements and for three Skyrme param- 
eterizations as indicated. 

hope for a similarly simple direct connection between fission barrier and fission life- 
time. A comparison of figure [3] with @] shows that the very gross trends may look 
similar. We have checked whether one could establish a simple relation between 
barriers and lifetimes. This turned out to be impossible. Fission life-times are very 
subtle quantities which depend on more ingredients than just barriers alone. 

There exist a few fission life-times along the decay channel from the upper island 
of stability. Fig. [6] compares computed life-times with data for fission and a-decay. 
The a life-times are well reproduced by the calculations, but the fission life-times 
are dramatically underestimated towards the lower end of the chain. The reason for 
that is still unclear. The mismatch calls for further investigations. 

4. Conclusions 

Fission life-times have been computed using the Skyrme-Hartree-Fock method. The 
fission path has been generated with a quadrupolc constraint where only axially 
symmetric deformations are considered. For the corresponding collective mass, we 
use self-consistent cranking. The quantum corrections to the collective potential 
(angular momentum projection, vibrational zero-point energy) are properly taken 
into account. The fission life-time is computed for thus given potential and mass by 
the WKB approximation, while the ground state energy, which is at the same time 
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Fig. 5. a-decay lifetime for a broad variety of super-heavy elements and for three Skyrme param- 
cterizations as indicated. 
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Fig. 6. Fission life-times and 
a-dccay life-times for 
elements along the decay 
chain from Z/N = 118/176 
computed with three dif- 
ferent Skyrme parameteri- 
zations, as indicated, and 
compared wit h experimental 
data fromEQl. 



the entrance energy for fission, is computed quantum mechanically. Results have 
been produced for three different Skyrme forces, SkP, SLy6 and SkI3, to explore 
the sensitivity to the parameterization. For comparison, we have also computed the 
a-decay life-times using the Viola systematics. 

A first test was performed by comparing with known fission life-times in the 
lower region of super-heavy elements, Z=104-108. The theoretical results gather 
around the experimental values with deviations of ±2 orders of magnitude. That 
can be called a promising agreement, the more so, as the isotropic trends are well 
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reproduced. The second test went for fission and a-decay life-times in the a-decay 
chain from 118/176. The a-decay life-times are well reproduced. But it turned out 
that the fission life-times are grossly underestimated by all three forces towards the 
lower end of the chain at 112/170. The reasons for that mismatch are not yet clear. 
It is to be reminded that these super-heavy nuclei have extremely soft deformation 
energy surfaces with several shape isomers which complicates the determination of 
the ground state. Notwithstanding that open problem, we have produced a sys- 
tematic survey of fission barriers, fission life-times and a-decay life-times for all 
conceivable super-heavy elements in the range Z=104-124 and N=150-190. Barriers 
and fission life-times show nicely the islands of shell stabilization around 104/152 
and 120/184. All forces predict a band of fission instability being orthogonal to the 
line connecting the two islands and crossing that around 112/166. There seems to 
be no path which could connect the 120/184 region with the lower island while 
avoiding fission. While the fission life-times show dramatic variation over the chart 
of super-heavy elements (from instability to 10 16 s), the a-decay times vary gently 
with small overall changes and without visible shell effects. The general crossover 
from a-decay to fission along the decay chains from the upper island is qualitatively 
reproduced by all three forces in the survey. 

Altogether, the results are promising and challenging at the same time. They 
call for further investigations exploring more systematically the sensitivity to the 
Skyrme parameterization and improving the description in the very soft transitional 
region around 112/166. 
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